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Abstract 

The target's characteristics play an important role in the 
design and performance analysis of all radar systems. The 
fluctuation rate of such a target may vary from essentially 
independent return amplitudes from pulse-to-pulse to 
significant variation only on a scan-to-scan basis. The 
correlation coefficient between the two consecutive echoes in 
the dwell-time is equal to unity for the SWI and SWIII cases 
while it is zero for the SWII and SWIV models. An important 
class of targets is that represented by the so-called 
moderately fluctuating Rayleigh and % 2 targets. The detection 
of partially correlated % 2 targets with two and four degrees of 
freedom is therefore of great importance. 

Our goal in this paper is to analyze the detection 
performance of the most familiar candidates of CFAR 
techniques for the case where the radar receiver post- 
detection integrates M pulses of an exponentially correlated 
signal from targets which exhibit % 2 statistics with two and 
four degrees of freedom. Exact formulas for the detection 
probabilities are derived, in the absence as well as in the 
presence of spurious targets. As predicted, the CA detector 
has the best homogeneous performance while the OS scheme 
gives the best target multiplicity performance when the 
number of outlying targets is within its allowable values. 
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I. Introduction 

The complete radar signal over a multiple observation 
interval is called a pulse train. Fluctuating pulse trains 
occur often in practice. When a radar target consists of 
several relatively strong reflecting surfaces displaced 
from one another by the order of a wavelength, the 
amplitude and phase of the composite radar echo are 
sensitive to the spatial orientation of the target. 
Moreover, if the target has relative motion with respect 
to the radar, it presents a time varying radar cross 



section. This change may be a slow variation and occur 
on a scan-to-scan basis (on successive antenna scans 
across a target) or on a pulse-to-pulse basis (on 
successive sweeps). Because it is difficult to predict the 
exact nature of the change, a statistical description is 
often adopted to characterize the target radar cross 
section [Meyer, 1973]. 

Target fluctuation is defined as variation in the 
amplitude of a target signal, caused by changes in 
target aspect angle, rotation, vibration of target 
scattering sources, or changes in radar wavelength. 
The fluctuation rate of a radar target may vary from 
essentially independent return amplitudes from pulse- 
to-pulse to significant variation only on a scan-to-scan 
basis. Three families of radar cross-section fluctuation 
models have been used to characterize most target 
populations of interest: the % 2 family, the Rice family, 
and the log-normal family [Aloisio, 1994]. The % 2 family 
includes the Rayleigh (SWI & SWII) model with two 
degree of freedom (k=1), the four degree of freedom 
(k=2) model (SWIII & SWIV), the Weinstock model 
(k<1), and the generalized model (k a positive real 
number). The y} models are used to represent complex 
targets such as aircraft and have the characteristic that 
the distribution is more concentrated about the mean 
as the value of the parameter k is increased [Swerling, 
1997]. The advantages of the Swerling target models are 
that they bracket a large number of real target classes 
[Di Vito, 1999]. 

It is often assumed that the Swerling cases bracket the 
behavior of fluctuating targets of practical interest. 
However, recent investigations of target cross section 
fluctuation statistics indicate that some targets may 
have probability of detection curves which lie 
considerably outside the range of cases which are 
satisfactorily bracketed by the Swerling cases. An 
important class of targets is represented by the so- 
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called moderately fluctuating Rayleigh and % 2 targets 
[El_Mashade, 2011] which when illuminated by a 
coherent pulse train, return a train of correlated pulses 
with a correlation coefficient in the range 0<p<l. The 
detection of this type of fluctuating targets is therefore 
of great interest [Hou, 1987]. 

The constant false alarm rate (CFAR) detection 
techniques represent the backbone upon which 
modern radar systems investigate in their operation. 
These processors make use of the fact that the 
amplitude variation of weather and sea clutter has a 
Rayleigh distribution, as well as the capability to 
reduce the clutter output to about the same level as the 
receiver noise level. As a consequence, much attention 
has been paid to the task of designing and assessing 
these adaptive detection techniques. The CA- and OS- 
CFAR detectors are the most widely used ones in the 
CFAR world [[El_Mashade, 2008]. While the 
performance of the CA detector is optimum in 
homogeneous situations, this performance degrades 
rapidly in non-ideal conditions caused by multiple 
target and non-uniform clutter. The OS processor, on 
the other hand, a CFAR technique, is relatively 
immune to non-homogeneous cases resulted from 
outlying targets and clutter edges. This technique relies 
on the order of the samples in the reference window 
and takes an appropriate reference cell to estimate the 
clutter power level. The OS trades a small loss in 
detection performance, relative to the CA, in ideal 
conditions for mush less performance degradation in 
non-ideal conditions. 

This paper is devoted to the analysis on the 
performance of CA- and OS-CFAR detectors for 
partially-correlated % 2 targets with two and four 
degrees of freedom in the absence as well as in the 
presence of spurious targets. In section II, we formulate 
the problem and compute the characteristic function of 
the post-detection integrator output for the case where 
the signal fluctuation obeys % 2 statistics with two and 
four degrees of freedom. The ideal performance of the 
schemes under consideration is analyzed in section III. 
Section IV deals with the problem of multiple-target 
environment and its effect on the behavior of the CA 
and OS processors. In section V, a brief discussion is 
made along with our conclusions. 

II. Theoretical Background and Statistical 
Model 

Let the input target signal and noise to the square-law 
detector be represented by the complex vectors u + jv 



and a + jb, respectively, u and v represent the in-phase 
and quadrature components of the target signal at the 
square-law detector, while a and b represent the in- 
phase and quadrature components of the noise, 
respectively. The target is assumed to be independent 
of the noise and the in-phase samples to be 
independent and identically distributed (IID) with the 
Gaussian probability density function (PDF), while the 
target samples are supposed to be identically 
distributed but correlated. The output of M-pulse non- 
coherent integrator, normalized to the noise power, is 
1 

\Ue ■ vk\ ■ \v i ' ut\ 



v A 



i ¥ S^ + aJ'+lv. + forj (1) 

Let the in-phase received target signal and noise 
vectors be U and A, respectively, where 

U A [a 1 ,a 2 ,a 3 , > u mY an< ^ 

A A [a 15 a 2 ,a 3 , ,a M J 

The characteristic function (CF) of v can be expressed 
as 

CS<°) ^ \p y (y) exp(-oy) dy 
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(3)In the above expression, pu (u) and pA (a) denote 
the joint probability density functions (PDF's) of U 
and A, respectively. Based on our previous 
assumptions, one can write the joint PDF of A as 



P A (a) 



exp 
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The substitution of Eq.(4) in Eq.(3) yields 
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There are many probability density functions for target 
cross section which are used to characterize fluctuating 
targets. The more important PDF is the so-called % 2 
distribution with 2k degrees of freedom. The % 2 family 
includes the Rayleigh (SWI & SWII) model, the four- 
degree of freedom model (SWIII & SWIV), the 
Weinstock model (k <1) and the generalized model (k a 
positive real number). These models used to represent 
complex targets such as aircraft have the characteristic 
that the distribution is more concentrated about the 
mean as the value of the parameter K is increased. The 
X 2 - distribution with 2k degrees of freedom is given by 



18 



Parallel and Cloud Computing Research (PCCR) Volume 1 Issue 2, July 2013 



www.seipub.org/pccr 



~o stands for the average cross section over all target 
fluctuations and U(.) denotes unit step function. When 
k=1, the PDF of Eq.(6) reduces to the exponential or 
Rayleigh power distribution that is applied to the 
Swerling case I. Swerling cases II, III, and IV 
correspond to k=M, 2, and 2M, respectively. When k 
tends to infinity, the j}' distribution corresponds to the 
non-fluctuating target. It is of importance to note that 
the % 2 - distribution with 2k degrees of freedom can be 
obtained by adding squared magnitude of k complex 
Gaussian random variables. 

II-l x 2 -Distribtion with Two-degrees of Freedom 

If k=1, then o may be generated as a=wi 2 +W2 2 , where 
Wi's are independent and identically distributed (IID) 
Gaussian random variables, each with zero mean and 
a/2 variance. The magnitude of the in-phase 
component u (u= I wi | ) has a PDF given by 

P (M) dtl = p (W, ) d W , = ^ 2y _ (7) 



exp 



2«r 



dw. with fix A 



cr 

~2 



To accommodate an Mxl vector of correlated y} RV's 
with two degrees of freedom, the PDF of the M- 
dimensional vector Wi is inreoduced. Therefore, the 

joint PDF of Ui's, i=l. 2, , M, has a form given by 

Prj^u ) du = P m (w,) dwx 



(8) 
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A is the correlation matrix of ui, U2, um and T 
denotes transposition. As I denotes the identity matrix, 
the substitution of Eq.(8) into Eq.(5) yields 

C^ = {^J * (9) 



If the average signal-to-noise ratio (SNR) is assumed to 
be a/ijj=Q, Eq.(9) can be re-expressed as 

c («) = n — 7—^ — (10) 

Ai's are the nonnegative eigenvalues of A. The CF of 
Eq.(lO) for the SWI target fluctuation model is 
represented by choosing Xi=M, Xi's=0, i=2, 3, . . ., M, and 
that for the SWII case can be modeled by the 
assumption that h's=l, i=l, 2, M. For the CFAR 
processor performance analysis to be simple, the above 
equation can be put in another clarified form as 



cm) =n 



J=1 co + d } 



& d, A 



i + ni 



(11) 



11-2 % 2 -Distribtion with Four-degrees of Freedom 



If K=2, then o and the squared magnitude of its in- 
phase component, u, may be generated as follows: 
Let w/s, i=l, 4, be IID Gaussian random variables, 
each one is of zero mean and of a= a/4 variance, a 
and | u | are given by [Weiner, 1988]: 



t 2 



and \u\ A ^w { 2 +w 2 



Since wi and W2 are IID, one may write 

p^fu) du = p (w,) p (w 2 ) dw x dw 2 
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Therefore, the joint PDF of ui, U2, , um is of the form 

P TT (u ) du = P (w,) (w 2 > dwidwi 



{— T 



exp 



A" 



W, + W 2 
2a 



A- 



(14) 



dw, dw 2 

where A is the correlation matrix of wi and W2. Note 
that wi and W2 are uncorrelated, but components of wi 
are correlated with each other, as the components of W2. 
Therefore, Eq.(5) can be rewritten as 

l 
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If the average SNR (2aA|/) is set as O, the above 
equation can be re-expressed as 

m + 1 (16) 



cm = n 



i + q4" \° j - 1 



The CF of Eq.(16) for the Swerling III target fluctuation 
model is represented by choosing Xi=M, Xis=0, i=2, 
3, M, and that for the Swerling IV case can be 
modeled by the assumption of ^i's=l, i=l, 2, M. To 
simplify the detection performance evaluation of the 
CFAR processors, the above equation can be 
reformatted to have an easily mathematical form such 
as 



(17) 



with b : A 



1 "H ClJl; /2 



In view of Eqs.(ll & 17), the solutionto partially- 
correlated case requires the computation of the eigen- 
values of the correlation matrix A. It is assumed here 
that i) the statistics of the signal are stationary, and ii) 
the signal can be represented by a first order Markov 
process. Under these assumptions, A is a Toeplitz 
nonnegative definite matrix. Thus, 
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Eqs.(ll, 17 &18) are the basic formulas of our analysis 
in this manuscript. 

The PDF of the output of the ith test tap is given by the 
Laplace inverse of Eq.(ll), in the case of % 2 distribution 
with two degrees of freedom, or Eq.(17), in the case of 
X 2 fluctuation with four degrees of freedom, after some 
minor modifications have been made. If the ith test tap 
contains noise alone,O=0 is set that is the average noise 
power at the receiver input \|/. If the ith range cell 
contains a return from the primary target, it rests 
without any modifications, where O represents the 
strength of the target return at the receiver input. On 
the other hand, if the ith test cell is corrupted by 
interfering target return, O must be replaced by I, 
where I denotes the interference-to-noise (INR) at the 
receiver input. 

III. Performance of CFAR Detectors in 
Benign Environments 

A benign environment in the present context is defined 
as an environment of additive Gaussian noise and 
possibly additive broadband jamming that can also be 
modeled as Gaussian. The statistical model with 
uniform clutter background describes the classical 
signal situation with stationary noise in the reference 
window. In this model, two signal cases are of interest: 
(a) one target in the test cell in front of an otherwise 
uniform background, (b) uniform noise situation 
throughout the reference window. In both cases, the 
noise neighborhood has a uniform statistic, i.e., the 

random variables Qi, Q2, , Qn in the reference 

window are assumed to be statistically independent 
and identically distributed. In the absence of the target 
case, the random variable v of the cell under test is 
assumed to be statistically independent of the 
neighborhood and subject to the same distribution as 
the random variables Qi's. 

The essence of CFAR is to compare the decision 
statistic v with an adaptive threshold TZ, see Fig. (1). 
the threshold coefficient T is a constant scale factor 
used to achieve a desired false alarm rate for a given 
window size N when the background noise is 



homogeneous. The statistic Z is a random variable 
whose distribution depends upon the particular CFAR 
scheme chosen and the underlying distribution of each 
of the reference range samples. Since the unknown 
noise power level estimate Z is a random variable, the 
processor performance is determined by calculating 
the average values of false alarm and detection 
probabilities. The false alarm probability Pfa is defined 
as 

P — ]p f M r tz) /v x. ,/x (19a) 

O 

with 

GO 

P fa (thr = a) A jp (x|Q = o) dx 



(19b) 



In the above expression, thr denotes the decision 
threshold and pz(.) represents the PDF of the noise 
power level estimate Z. In terms of the CF of the 
random variable Z, its PDF can be obtained by taking 
the inverse Laplace transformation of this CF. Thus, 

n (y) = t^ICz^) exp(«y) deo (20) 

= 22 res { C z C<«) exp(ai y), co j } 

j 

The contour of integration 'c' lies in the right of all 
singularities of Cz(ro) beginning in the left half-plane, 
coj's are the poles of Cz(a>), and res(.) denotes the 
residue. Hence, 



f p (.y) dy = -J^res \ < -- A 1 , r>, 



(21) 



Finally, the substitution of Eqs.(21, 20 & 19) into the 
definition of Pfa results in 



CO 



(22) 



Similarly, the detection probability can be easily 
calculated as 

v f CvH"*°) C z (-^> I (23) 

The detection probability of a CFAR processor for % 2 
targets with two degrees of freedom can be obtained 
by substituting Eq.(ll) into the definition of Pd, which 
gives 

P = - > res{ — ■ — , <y 



(24) 



i n 



On the other hand, when the fluctuation of the primary 
target follows the x 2 distribution with four degrees of 
freedom, C v (ro) takes the following simplified form 

<° + 1 (25) 



where 
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The substitution of Eq.(25) into the definition of the 
detection probability gives the processor detection 
performance for % 2 targets with four degrees of 
freedom, which becomes 



(28) 



To obtain the relationship between the false alarm 
probability Pfa and the thresholding constant T, it is 
assumed that no target is present in the test cell (A=0). 
Therefore, as M-pulse non-coherent integration is used, 
Pfa can be calculated as 

P f .-^^\C^\ .(29, 



r(Af) daf 



CO 



co=-\ 



In this case, C v (ro) has an Mth order pole at &=-\. Thus, 

^"§ix7M)^ ( - r -V-i (30> 

It is of important to note that the above expression for 
the rate of false alarm is verified either for % 2 targets 
with two degrees of freedom or for x 2 targets with four 
degrees of freedom. Moreover, the characteristic 
function of the noise power level estimate 'Z' is the 
backbone of the processor performance analysis, as 
shown in Eqs.(24. 28 & 30), either in homogeneous or 
non-homogeneous background environments. 
Therefore, our scope in the following subsections is to 
calculate this important quantity for the two processors 
under consideration. 

a) Cell-Averaging (CA) Detector 

A simple approach to achieve the CFAR condition is to 
set the detection threshold on the basis of the average 
noise power in a given number of reference cells where 
each of these cells is assumed to contain no targets. 
Such a scheme is denoted as cell-averaging (CA) CFAR 
processor. This detector is specifically tailored to 
provide good estimates of the noise power in the 
exponential PDF. In this CFAR detection technique, the 
total noise power level is estimated by the sum of N 
range cells of the reference window. This is a complete, 
sufficient statistic for the noise power \|/ under the 
assumption of exponentially distributed homogeneous 
background. In other words, the noise power is 
estimated as 



and 



(31) 



Q, 



Z a AEI q u aZ Q, 

i=i j=i 1 i=i 

M 

Under the assumption that the surrounding range cells 
contain independent Gaussian noise samples with the 
same variance, Zca is the maximum likelihood estimate 
of the common variance. In the absence of the radar 
target, each cell of the reference window has a CF of 
the same form as that of v with the exception that O 
must be set to be zero. Therefore, the elimination of O 
in either Eq.(ll) or Eq.(17) leads to 

1 



CO 



+ 1 



(32) 



Since the random variables Qi's are assumed to be 
statistically independent, ZCA has a CF given by 

f 1 \ L 



C 7 (®) = 
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CO + l 
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(33) 



Once the CF of the noise power level estimate is 
calculated, the detector performance evaluation can be 
easily carried out when the primary target fluctuates in 
accordance with j2 model either with two or four 
degrees of freedom as shown in Eqs.(24, 28 & 30), 
where 

| ( i \ L+K 

(aa))\=(-a) K (L) K 



da 



aco +1 



J 



(34)The Pochammer symbol (x)X is as defined in 
[[El_Mashade, 1994]]: 

(X) A (35) 
y J> ~ \x{x+\) (x+l-l) for £=1,2,3,... 

The rationale for the CA type of CFAR schemes is that 
by choosing the mean, the optimum CFAR processor 
in a homogeneous background when the reference 
cells contain IID observations governed by exponential 
distribution is achieved. As the size of the reference 
window increases, the detection probability 
approaches that of the optimum detector which is fixed 
threshold architecture [[El_Mashade, 2006]. 

b) Order-Statistics (OS) Detector 

The performance of CA-CFAR detector degrades 
rapidly in non-ideal conditions caused by multiple 
targets and non-uniform clutter. The ordered-statistic 
(OS) CFAR is an alternative to the CA processor. The 
OS trades a small loss in detection performance, 
relative to the CA scheme, in ideal conditions for much 
less performance degradation in non-homogeneous 
background environments. 
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Order statistics characterize amplitude information by 
ranking observations in which differently ranked 
outputs can estimate different statistical properties of 
the distribution from which they stem. The order 
statistic corresponding to a rank K is found by taking 

the set of N observations Qi, Q2 , Qn and ordering 

them with respect to increasing magnitude in such a 
way that 

e (1) <-<e ( „)<e m <e (M) <-<e w 

where Q(K) is the Kth order statistic. The central idea of 
OS-CFAR procedure is to select one certain value from 
the above sequence and then used as an estimate Z for 
the average clutter power as observed in the reference 
window. Thus, 



ws " i-'-'-' >-■ ) (37) 

The OS scheme with parameter K will be denoted as 
OS(K). The value of K is generally chosen in such a 
way that the detection of radar target in homogeneous 
background environment is maximized. 
In order to analyze the processor detection 
performance in uniform clutter background, the CF of 
the random variable Zos is required. This CF can be 
expressed as [[El_Mashade, 2011]] 

K ( N %A ( K- 
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(40)Once the CF of the noise power level estimate Zos 
is computed, the processor detection performance 
becomes an easy task, as we have previously proved. 
Finally, the €th derivative of this CF is given by: 

(41) 
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A desirable CFAR scheme would of course be one that 
is insensitive to changes in the total noise power within 



the reference window cells so that a constant false 
alarm rate is maintained. This is actually the case of the 
two architectures under consideration. 

c) Numerical Results 

In this section, we present some representative 
numerical results, which will give an indication of the 
tightness of our previous analytical expressions. The 
performance of CFAR processors for partially- 
correlated x 2 fluctuating targets is numerically 
evaluated for some parameter values and the results of 
these evaluations are presented in several sets of 
figures. The first set includes Figs.(2-3) and concerns 
with the detection performance of CA and OS schemes, 
respectively, for a number of integrated pulses of 2, 3, 
and 4 along with the single-sweep case which is 
included as a reference for comparison, when the radar 
target fluctuates in accordance with partially- 
correlated x 2 model with two and four degrees of 
freedom. Since the performance of OS processor is 
strongly dependent on the ranking order parameter K, 
we choose the value that corresponds to the optimum 
detection performance in uniform noise background 
which is 21 for N=24 [8]. The displayed results of these 
figures show that for low values of SNR, the fully- 
correlated case (# 1 =1.0) gives higher detection 
performance than the fully de-correlated (# 1 =0) case. 
As the target return becomes stronger, an alternative 
version of the above behavior is observed, where the 
non-correlated performance surpasses the fully- 
correlated one. It is also noted that for M=2, the 
processor detection performance for ® i & ® 1 =0 
coincides with its detection performance for ® 3 & ® 1 
=1.0. In addition, the non-correlated and fully- 
correlated detection performances embrace all the 
partially-correlated cases in either situation. This 
observation is common for any number of integrated 
pulses. Additionally, the partially-correlated % 2 with 
© 3 offers higher detection performance than that 
with # 4 when the target return becomes stronger 
while the reverse of this behavior occurs when the 
target return is modest. Moreover, as M increases, the 
processor detection performance ameliorates and the 
reversing point is shifted towards lower values of SNR. 
By reversing point, we mean the point at which the 
detection performance changes its superiority from 
fully-correlated to fully-uncorrelated. All the presented 
results are obtained for a constant false alarm rate of 
10 6 and a reference window of size 24 cells. 

To make a comparison between the performances of 
the underlined CFAR schemes and that of the 
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optimum detector, the second group of illustrations 
includes Figs. (4- 7). This category contains the partially- 
correlated x 2 target homogeneous detection 
performance of the CA and OS schemes along with the 
optimum processor for M=2 and 4. The curves of this 
set are labeled in the CFAR procedure and the 
correlation coefficient 'qi', respectively. The indication 
OT, A, or O on a specified curve means that it is drawn 
for optimum, cell-averaging, or order-statistic detector, 
respectively. It is of importance to note that Figs.(4-5) 
describe the performance of the three processors when 
the primary target fluctuates following partially- 
correlated x 2 model with two and four degrees of 
freedom, respectively, for M=2, while Figs. (6-7) depict 
the same thing for M=4. In any situation, it is noticed 
that the performance of CA algorithm is the much 
closer one to that of the optimum detector under any 
operating conditions and the OS scheme comes next. 
All the indicated remarks on the results of Figs.(2-3) 
are clearly demonstrated on the results of the present 
figures. For small SNR, the processor performance 
degrades as qi decreases, while for large SNR, this 
performance improves as qi decreases. This 
observation is noticed for partially-correlated x 2 targets 
fluctuating with either two or four degrees of freedom. 
In addition, these plots illustrate the superiority of x 2 
fluctuation model with 4 degrees of freedom over that 
following % 2 model with 2 degrees of freedom, 
especially for large SNR. As the number of non- 
coherently integrated pulses increases, the processor 
performance improves and less SNR value is needed to 
achieve the same level of detection. 

In the next group of figures, we are concerned with 
what is known, in the world of radar target detection, 
as receiver operating characteristics (ROC's). These 
characteristics describe the detection probability as a 
function of the false alarm probability for a fixed target 
signal strength (SNR=const). This set of figures 
incorporates Figs. (8-9) which represent the ROC's of 
CA and OS(21), respectively. In these figures, the 
processor detection performance is plotted against the 
false alarm rate for a number of integrated pulses of 2 
and 4 when the primary target fluctuates in accordance 
with y} distribution with two and four degrees of 
freedom. Since the fully correlated and de-correlated 
cases enclose the partially-correlated situations, our 
evaluations are restricted to these two border limits 
(m i=0 & # i=l). The reference window size is taken as 
24 and the primary target signal strength as 5 dB. For 
comparison, the single sweep ROC's is also included in 



these plots. The curves of these figures are labeled in 
the number of post-detection integrated pulses and the 
correlation coefficient of the primary target returns 
(# i). The displayed results show that for lower false 
alarm rates, the processor detection probability 
improves as # i increases. As the false alarm rate 
increases, this behavior is rapidly reversed and the 
fully de-correlated detection performance surpasses 
that corresponding to the fully correlated case. 
Additionally, as the number of non-coherently 
integrated pulses increases, the reversing point or the 
critical rate, moves towards the lower false alarm rate 
and this is common as the target fluctuates following % 2 
with either two or four degrees of freedom. Moreover, 
the processor detection performance for % 2 distribution 
with ® =2 is supeior to that for x 2 distribution with 
# =1 in the case where the operating false alarm rate is 
higher than its critical value, otherwise, the processor 
detection performance when the fluctuation of the 
primary target follows x 2 distribution with ® =1 
surpasses its performance when this fluctuation obeys 
X 2 distribution with # =2Keeping in mind that when 
M=2, the processor ROC's for ® 4 & ® i=0 coincide 
with its ROC's for ® 2 & © i=l. The results of this 
group confirm the observed remarks on the previous 
groups as well as demonstrate our comments on the 
behavior of their figures. Additionally, the 
homogeneous detection performance of CA scheme is 
always superior to that of OS algorithm, under the 
same operating conditions, as predicted. 

To illustrate the influence of the signal correlation on 
the processor detection performance, the last group is 
devoted to the processor detection performance 
against the correlation coefficient between the radar 
target returns given that the number of integrated 
pulses (M) as well as the signal strength (SNR) is held 
constant. Figs.(lO-ll) illustrate these characteristics for 
the underlined detectors, CA and OS(21), respectively, 
for SNR=-5, 0, and 5dB when M=l, 2, and 4, given that 
the radar target fluctuates in accordance with % 2 
distribution with <s> =1 Since the chosen signal strength 
is modest, the scheme detection performance improves 
as the target returns become highly correlated. This 
result is predicted since the correlation strengthens the 
weak signal returns in contrast to the case in which the 
signal returns are strengthened but the correlation 
weakens the signal returns. For the same reason, the 
rate of improvement decreases as the number of 
consecutive sweeps increases. Moreover, this rate of 
improvement increases as the target signal returns 
become weaker, given that the number of integrated 
pulses maintains constant. The two processors under 
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consideration give approximately the same level of 
detection for correlated signal returns without any 
superiority between them. This concluded remark is 
also predicted due to the fact that the CFAR schemes 
behave likely for the modest target signal strength, as 
demonstrated by the above mentioned results. Finally, 
if the primary target fluctuates in accordance with % 2 
distribution with ® =2,the underlined processors act 
the same behavior with minor improvement in each 
case. 

IV. Performance of CFAR Schemes in 
Multitarget Environment 

The performance of the CFAR algorithms for uniform 
clutter model is completely evaluated in the previous 
section. Clutter edges, on the other hand, are used to 
describe transition areas between regions with very 
different noise characteristics. Since it is concerned 
here with partially correlated x 2 targets, this situation 
is of secondary scope. On the other hand, multiple 
target situations occur occasionally in radar signal 
processing when two or more targets are at a very 
similar range. The consequent masking of one target by 
the other is called suppression. These interferers can 
arise from either real object returns or pulsed noise 
jamming. From a statistical point of view, this implies 
that the reference samples, although still independent 
of one another, are no longer identically distributed. 
Let us now examine the dependence of the 
performance of the CFAR procedures on the accurate 
knowledge of the target fluctuation model when the 
reference window is contaminated with a fluctuating 
interfering target returns. In our study of the non- 
homogeneous background, the amplitudes of all the 
targets present amongst the candidates of the reference 
window are assumed to be of the same strength and to 
fluctuate in accordance with the partially-correlated y} 
fluctuation model with correlation coefficient Qi. The 
interference-to-noise ratio (INR) for each of the 
spurious targets taken as a common parameter is 
denoted by I. 

a) Cell-Averaging (CA) Detector 

To analyze the CA performance when the reference 
window no longer contains radar returns from a 
homogeneous background, the assumption of 
statistical independence of the reference cells is 
retained. Suppose that the reference window contains r 
cells from interfering target returns with background 
power of ^(1+1) and N-r cells from clear background 



with noise power \\>. Thus, the estimated total noise 
power level is obtained from 

z CA = t ^ + S q, — x + y < 43 ) 

The random variable representing the interfering 
target return 'X' has a CF of the same form as that 
given by Eq.(ll), when the extraneous targets fluctuate 
in accordance with x 2 model with two-degrees of 
freedom, and as that given by Eq.(16), in the case 
where the spurious targets are of the % 2 model with 
four-degrees of freedom. In both cases, A (SNR of the 
primary target) should be replaced by I (INR of the 
secondary target). On the other hand, since each RV in 
Y represents the clear background return, it has a 
similar form, for its CF, as given by Eq.(32). Since the 
candidates of each type are assumed to be statistically 
independent, we have 

c z J°» = \n^tr\ \—A & (44) 
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In the above expression, Ai's represent the nonnegative 
eigenvalues of the correlation matrix, with Q=Qi in Eq. 
(18), associated with the spurious target returns. 
Finally, the derivative of the CF of the test statistic of 
the CA procedure, when the interfering targets 
fluctuate in accordance with % 2 model of four degrees 
of freedom, can be easily computed as 



(46) 
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We repeat again that once CF of the noise power level, 
Zca, is obtained, the processor performance evaluation 
is completely determined, in the case where the 
fluctuation of the spurious targets follows the % 2 model 
with either two or four degrees of freedom, as 
previously demonstrated. 

b) Ordered-Statistics (OS) Detector 

In order to analyze the processor performance when 
there are interfering target returns amongst the 
contents of the reference window, the assumption of 
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statistical independence of the reference cells is 
retained. Consider the situation where there are r 
reference samples contaminated by extraneous target 
returns, each with power level ^(1+1), and the 
remaining N-r reference cells contain thermal noise 
only with power level t|;. Under these assumptions, the 
Kth ordered sample, which represents the noise power 
level estimate in the OS detector, has a cumulative 
distribution function (CDF) given by [El_Mashade, 
2001]: 



N mm(i,N-r) 

I I 

i=K j=max(0,/-j-) 



(N-A 






v-jj 



F z (*) = ! I 

£-> OS 

b-FSzt r - j [i-FM r - i+ii {F s (zt j 



(47) 



The CDF of the reference cell that contains a spurious 
target return, when this target fluctuates in accordance 
with x 2 with two degrees of freedom, can be calculated 
as: 
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On the other hand, if the interfering target's fluctuation 
follows chi-square model with four degrees of freedom, 
Fs(z) takes the form 
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Since F(x)=l-[1-F(x)], Eq.(47) can be written in another 
simpler form as [15]: 
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By taking the Laplace inverse of Eq.(32) and then 
substitued, along with Eq.(48), into Eq.(53) one obtains, 
for % 2 target fluctuation with two degrees of freedom, 
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On the other hand, if the extraneous targets fluctuate 
according to % 2 model with four degrees of freedom, 
Eq.(53) takes the form 
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To determine the detection performance of the OS- 
CFAR processor, it is important to calculate the 
Laplace transform for its test statistic, where the false 
alarm and detection probabilities are completely 
dependent on this transformation along with its 
derivatives with respect to co. The co-domain 
representation of Eq.(54) is: 
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On the other hand, the Laplace transformation of 
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Eq.(55) is analytically very complicated and therefore 
the numerical techniques is the only way that allows us 
to compute this transformation. 

Again, the OS-CFAR processor performance is highly 
dependent upon the value of K. For example, if a single 
extraneous target appears in the reference window of 
appreciable magnitude, it occupies the highest ranked 
cell with high probability. If K is chosen to be N, the 
estimate will almost always set the threshold based on 
the value of interfering target. This increases the overall 
threshold and may lead to a target miss. If, on the other 
hand, K is chosen to be less than the maximum value, 
the OS-CFAR scheme will be influenced only slightly for 
up to N-K spurious targets. 

c) Numerical Results 

Here, we provide a variety of numerical results for the 
performance of CA- and OS-CFAR processors in 
multiple target environments. Since the optimum value 
of K, for a reference window of size 24, is 21, it has been 
assumed that there are three interfering target returns 
amongst the contents of the estimation cells. This value 
of extraneous target returns is the maximum allowable 
value before the OS performance degradation occurs. 
Our numerical results are obtained for a possible 
practical situation where the primary and the secondary 
interfering targets fluctuate in accordance with the % 2 
fluctuation model with the same correlation coefficient 
(qi=02=# ), and of equal target return strength 
(INR=SNR). Additionally, for our new results to be 
comparable with the older ones (in the absence of 
outlying targets), the design rate of false alarm is held 
unchanged (Pfa=10 6 ). The obtained results are classified 
into categories. The first category, which includes 
Figs.(12,13), depicts the detection performance of the CA 
and OS(21) schemes, respectively, under the same 
operating conditions as in Figs. (2,3) with the exception 
that three spurious targets are allowed to be present 
amongst the candidates of the reference set. The curves 
of these figures are labeled in the number of integrated 
pulses (M), the strength of correlation between the 
target returns (® )and the parameter that represents the 
degrees of freedom (# ).For comparison, the single 
sweep case is included in these plots. In addition to the 
partially-correlated cases, the well-known four 
Swerling's models are also included. It is important to 
note that the full scale of the Pd axis in this case of CA 
scheme has a maximum value of 75%. There is another 
interesting point about the reaction of CFAR schemes 
against % 2 fluctuation model where the SWII & SWIII 
performances are coincide for M=2. The displayed 



results show that the behavior of the CFAR processors 
in the presence of extraneous targets is similar to that in 
their absence. By comparing the results of Fig.(12) to the 
corresponding ones of Fig.(13), it is observed that 
intolerable masking of the primary target occurs in the 
case of CA procedure and the OS(21) architecture is 
capable to resolve multiple targets in the reference 
window as long as their number is within its allowable 
values (r < N-K). All the concluded remarks about the 
ideal performance of CA and OS(21) schemes are also 
observed on their multi-target performance given that 
the target under test along with the interfering ones 
fluctuate obeying x 2 model with 2 and 4 degrees of 
freedom. For the comparison between the behaviors of 
the considered detectors, against the presence of 
outlying targets to be clarified, the next category 
combines their performances in the same figure. This 
category, containing Figs.(14-15), illustrates the 
multiple- target detection performance of CA and OS(21) 
processors when the principal and the spurious targets 
fluctuate following x 2 model with 2 and 4 degrees of 
freedom, respectively. As a reference level of 
performance, the detection performance of the optimum 
scheme is also included in these figures. The superiority 
of the OS performance is very clear under any level of 
correlation and the x 2 fluctuation model with 4 degrees 
of freedom gives always an improved performance than 
that fluctuating with 2 degrees of freedom. This result is 
predicted since all the reference cells participate in 
constructing the detection threshold in the case of CA 
processor while the K th sample is the only one that is 
responsible for building the detection decision. This 
means that the presence of interferers will raise the 
threshold in the case of CA scheme while it has no effect 
on the detection threshold in the OS algorithm, given 
that their number is within its allowable values. 
Additionally, the displayed results confirm the 
improvement of the processor detection performance 
with increasing the number of post-detection integrated 
pulses, as the theoretical analysis demonstrates. The 
next set of curves (Figs.16-17) illustrates the required 
SNR to achieve an operating point (Pfa=10 6 j Pd=90%) for 
the two CFAR processors under investigation, as a 
function of the consecutive pulses correlation coefficient 
(01=02=0), for several values of post-detection integrated 
pulses when the fluctuation model of the tested target as 
well as the secondary interfering ones obeys % 2 
distribution with 2 and 4 degrees of freedom, 
respectively. The candidates of this set are labeled in the 
CFAR processor, number of integrated pulses, and the 
number of reference cells that are contaminated with 
extraneous target returns. It is of importance to note 
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that r=0 means homogeneous background which is free 
from any interferers. In that case, the CA processor 
needs less strengthened target signal compared to the 
OS scheme, to attain the required level of detection for 
a given rate of false alarm. In addition, the required 
SNR increases as q increases, given that M along with r 
are held fixed. In multi-target situation, on the other 
hand, the CA procedure is unable to reach the given 
level of detection (90%) for the specified values of its 
parameters. For this reason, Fig.(16) shows the 
multiple-target detection performance of the OS 
processor only since the CA scheme is unable to 
achieve a detection probability of 0.9 value in the 
presence of 3 outlying targets, for either M=2 or M=4. 
Since the required value of the probability of detection 
is high, the processor detection performance for % 2 
fluctuating targets with 4 degrees of freedom is better 
than its performance for those fluctuating in 
accordance with % 2 model with 2 degrees of freedom. 
This is the fact that Fig.(17) demonstrates. Finally, the 
variations of the false alarm rate with the strength of 
the interfering target returns are plotted, in the last 
category of curves, for the underlined architectures 
when they are reacted to % 2 fluctuating targets with 2 
and 4 degrees of freedom. This set of figures includes 
Figs. (18-19). In this case, the primary and the 
secondary outlying targets are assumed to be 
fluctuating with the same degree of freedom and our 
numerical results are constrained to the two limits of 
the correlation coefficient (o=0, Q=l) which are 
corresponding to the well-known Swerling's models. 
The curves of these figures are labeled in the specified 
processor (A: for cell- Averaging, O: for Order-statistic), 
number of integrated pulses M and the Swerling 
fluctuation model. The numerical results of these 
figures are given for a design rate of false alarm of lO 6 
and a reference window of size 24 samples. The 
interfering targets are assumed to exhibit y} statistics 
with full correlation (oi=l), which corresponds to SWI, 
when the degrees of freedom is 2, or SWIII, when the 
degrees of freedom is 4. On the other hand, when the 
outlying targets exhibit % 2 statistics with null 
correlation (Qi=0), the resulting fluctuation model is 
SWII, in the case where the degrees of freedom is 2, or 
SWIV, in the case where the degrees of freedom is 4. 
The results of these figures show that the false alarm 
rate performance of the CA processor degrades as the 
strength of the interfering target return (INR) increases, 
and the rate of degradation decreases as o increases or 
the degree of freedom decreases. This statement is 
common for the CFAR processors considered here. 
However, the OS(21) scheme has the lowest rate of 



degradation and consequently, it is the only processor 
capable to maintain the false alarm rate at 
approximately the desired value, given that the 
number of spurious target returns doesn't exceed its 
allowable values. In addition, when the INR becomes 
very high, the behavior of the OS detector is 
independent on the correlation coefficient q and tends 
to be constant. Moreover, the OS false alarm rate 
performance improves as the number of post-detection 
integrated pulses increases, which is excluded from 
offer in the CA technique. This result is expected since 
the largest interfering target returns occupy the top 
ranked cells and therefore they are not incorporated in 
the estimation of the background noise power level. In 
other words, the noise estimate is free of extraneous 
target returns and therefore it represents the 
homogeneous background environment. As a result, 
the false alarm rate performance of the OS procedure 
improves as the number of non- coherently integrated 
pulses increases. 

V. Conclusions 

In this paper, the detection probability of a radar 
receiver that post-detection integrates M pulses of an 
exponentially correlated signal from a Rayleigh target 
in thermal noise is calculated. At the limiting 
correlation coefficients, p=l and p=0, the analysis yields, 
respectively, the well-known Swerling cases 1 & 2. In 
addition, the probability of detection of the sum of M 
square-law detected pulses is computed for the case 
where the signal fluctuation obeys x 2 statistics with 
four degrees of freedom. Swerling's well-known cases 
3 & 4 represent the situations where the signal is 
completely correlated and de-correlated, respectively, 
from pulse to pulse. Moreover, the performance of 
CFAR processors for partial signal correlation in these 
two situations has been analyzed. These processors 
include CA detectorwith the best homogeneous 
performance, and OS procedure, which trades a small 
loss in homogeneous performance, relative to CA, for 
much less performance degradation in target 
multiplicity environments. The exact false alarm and 
detection probabilities have been derived, in the 
absence as well as in the presence of extraneous targets, 
for the condition of partial signal correlation. The 
primary and the secondary interfering targets are 
assumed to be fluctuating in accordance with % 2 - 
distribution with two and four degrees of freedom. The 
analytical results have been used to develop a 
complete set of performance curves including ROC's, 
detection probability in homogeneous and multi-target 
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situations, required SNR to achieve a prescribed values 
of Pfa & Pd, and the variation of false alarm rate with 
the strength of interfering target returns that may exist 
amongst the contents of the estimation set. As expected, 
the detection performance of the CFAR detectors for 
partially-correlated x 2 targets with two degrees of 
freedom is between that for SWI targets and that for 
SWII models. On the other hand, the detection 
performance of the CFAR schemes for partially 
correlated x 2 targets with four degrees of freedom lies 
between that for SWIII fluctuation model and that for 
SWIV targets. In any one of these fluctuating families, 
more per pulse signal-to-noise ratio is required to 
achieve a prescribed probability of detection as the 
signal correlation increases from zero to unity. In 
addition, the false alarm rate increases with the signal 
correlation and the OS architecture is the only, relative 
to the CA scheme, processor that is capable to maintain 
a constant rate of false alarm, irrespective to the 
interference level, in the case where the spurious target 
returns occupy the top ranked cells and they are within 
their allowable values. As a final conclusion, the 
detection performance of a CFAR procedure is related 
to the target model, the number of post-detection 
integrated pulses, and the average power of the target. 
When the target signal fluctuates obeying % 2 statistics, 
the signal components are correlated from pulse to 
pulse and this correlation degrades the processor 
performance. A common and accepted practice in 
radar system design to mitigate the effect of target 
fluctuation is to provide frequency diversity to de- 
correlate the signal from pulse to pulse. While this 
technique is effective, it requires additional system 
complexity and cost. 

The results for partial-correlation fall between those for 
the two extremes of complete de-correlation and 
complete correlation. Thus, to estimate the 
performance for partially-correlated pulses, 
interpolation between the results for completely 
correlated and de-correlated conditions can be used as 
an approximation. Extension to the partially-correlated 
%2 target with 2k degrees of freedom, where k>2, is 
relatively straightforward. 
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Fig.(2) M-sweeps ideal detection performance of CA scheme for partially-correlated %2 targets 
with 2 & 4 degrees of freedom when N=24, and Pfa=1.0E-6. 
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Fig. (3) M-sweeps ideal detection performance of OS(21) scheme for partially-correlated %2 targets 
with 2 & 4 degrees of freedom when N=24, and Pfa=1.0E-6. 
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Fig. (4) M-sweeps ideal detction performance of CFAR processors for partially-correlated %2 targets 
with two degrees of freedom when N=24, M=2, Pfa=1 .OE-6. 
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Fig.(5) M-sweeps ideal detction performance of CFAR processors for partially-correlated %2 targets 
with four degrees of freedom when N=24, M=2, Pfa=1.0E-6. 
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Fig.(6) M-sweeps ideal detction performance of CFAR processors for partially-correlated %2 
targets with two degrees of freedom when N=24, M=4, Pfa=1.0E-6. 
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Fig.(7) M-sweeps ideal detction performance of CFAR processors for partially-correlated %2 
targets with four degrees of freedom when N=24, M=4, Pfa=1 .OE-6. 
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Fig.(8) M-sweeps receiver operating characteristics (ROC's) of CA processor for partially-correlated 
%2 targets in homogeneous situation when N=24 and Q=5dB. 
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Fig.(9) M-sweeps ideal receiver operating characteristics (ROC's) of 0S(21) scheme for 
partially-correlated %2 targets when N=24, and fi=5dB. 
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Fig. (10) M-sweeps variation of ideal detection performance of CA-CFAR processor for partially- 
correlated Rayleigh fluctuating targets when N=24, Pfa=1.0E-6, and SNR=-5, 0, & 5dB. 
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Fig.(11) M-sweeps variation of ideal detection performance of OS(21)-CFAR processor for partially- 
correlated Rayleigh fluctuating targets when N=24, Pfa=1.0E-6, and SNR=-5, 0, & 5dB. 
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Fig.(12) Multitarget detection performance of CA scheme for partially-correlated %2 targets 
with 2 & 4 degrees of freedom when N=24, r=3, and Pfa=1 .OE-6. 
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Fig.(13) Multitarget detection performance of OS(21) scheme for partially-correlated %2 targets 
with 2 & 4 degrees of freedom when N=24, r=3, and Pfa=1.0E-6. 
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Fig.(14) M-sweeps multitarget detction performance of CFAR processors for partially-correlated %2 
targets with two degrees of freedom when N=24, M=3, r=3, Pfa=1.0E-6. 
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Fig.(15) M-sweeps multitarget detction performance of CFAR processors for partially-correlated %2 
targets with four degrees of freedom when N=24, M=3, r=3, Pfa=1.0E-6. 
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Correlation coefficient "p " 

Fig.(16) M-sweeps required SNR, in the absence as well as in the presence of outlying targets, 
to achieve an operating point (1.0E-6, 0.9) of CFAR schemes for partially-correlated %2 targets 

when N=24, and k=2. 
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Fig. (17) M-sweeps required SNR, in the absence as well as in the presence of outlying targets, to 
achieve an operating point (1.0E-6, 0.90), of CFAR schemes for partially-correlated %2 targets when 

N=24, and k=4. 
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Interference-to-noise ratio (INR) in "dB" 

Fig. (19) M-sweeps false alarm rate performance, as a function of the interfering target strength 
(INR), of CFAR schemes for %2 fluctuating targets when N=24, r=3, k=4, and design Pfa=1.0E-6. 



39 



